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Abstract 



Dissipative systems can be described in terms of non-hermitian hamil- 
tonians H, whose left eigenvectors and right eigenvectors |/,) form a 
bi-orthogonal system. Bi-orthogonal systems could suffer from two difficul- 
ties, (a) If the eigenvectors do not span the whole space, then H can only be 
diagonalized to blocks (the Jordan-block problem), (b) Normalization would 
not be possible and many familiar-looking formulas would fail if (f^lfj) = 
for some j (the orthonormalization problem). Waves in open systems provide 
a well-founded realization of a bi-orthogonal system, and it is shown that 
these two problems can indeed occur and are both related to higher-order 
poles in the frequency-domain Green's function. The resolution is then given 
by introducing a generalized duality transformation involving extra basis vec- 
tors, whose time evolution is modified by polynomials in the time t. One thus 
obtains a nontrivial extension of the bi-orthogonal formalism for dissipative 
systems. 

PACS numbers: 03.40.Kf, 02.30.Mv, 02.60.Lj, 03.65.-w 
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I. INTRODUCTION 



A. Dissipative systems and bi-orthogonal bases 

Dissipative systems are often discussed, in a phenomenological way, by postulating a non- 
hermitian hamiltonian (NHH) H, whose left eigenvectors and right eigenvectors \fj) form 
a bi-orthogonal system (BS) Jl|fj,0,||,0 . It is usually assumed that these eigenvectors are 
complete; the BS then constitutes a bi-orthogonal basis (BB). These NHHs with discrete 
BBs can sometimes be obtained from a full microscopic theory, but usually under some 
approximations [|3],^| . BBs allow dissipative systems to be placed into a familiar framework; 
the advantages are obvious and need not be enumerated. As a consequence, there is a 
substantial literature on both the mathematics and the physical applications of BBs (e.g., 
(chemical bonding), M, M (solid mechanics)). 

To introduce the formalism consider a Hilbert space with an NHH H, supposed to 
have a complete basis of eigenvectors \fj), with eigenvalues Uj. In terms of the standard 
inner product (-|-), introduce the unique dual basis {I/- 7 )} with 

(f j \fk) = S jk ; (1.1) 

one finds (f k \W\p} = (f j \H\f k )* = oo* k 6 jk , so that = oj*\f j ). Now define a duality 

transformation V by T>\fj) = \P), extended to the whole space by conjugate linearity [JTU| : 



V («|X) + m) = otV\ X ) + FVty) . (1.2) 

The symmetry of ( |1 . 1| ) shows that the duality transformation satisfies T> 2 = 1, justifying its 
name 0. Since (f k \VH\fj) = u*5 jk = (M^Vlfj), one has 

VH = E X V . (1.3) 
Any vector can now be expanded as 

10) = (i.4) 

j 

{P\4>) 

Jim- (L5) 

where for application in the following Sections it is convenient not to fix the normalization 
(f^lfj) by (|1.1|) . This leads immediately to the resolution of the identity and of the time- 
evolution operator 



r (/' ./}> 



e 



-iHt 



v l/i)e- Wjt </ J l ri7 x 



which in principle solves all the dynamics. For example, the wavefunction at time t > 
with (11.41) as initial data would be 
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|0(t)>=J>,e-^|/,> 
j 



(1.8) 



It is emphasized that (|1.3| )- (|1.8| ) have been obtained only for the case in which {\fj)} con- 
stitutes a complete basis. 

This brief account exposes two related lurking difficulties that are usually brushed aside: 
the Jordan block problem and the orthonormalization problem. 



B. Jordan blocks and orthonormalization 

Because H is not hermitian, there is no guarantee that its eigenvectors \fj) form a 
complete set, so a BB may not exist. There are in fact two fundamentally different reasons 
why the eigenvectors may fail to be complete [12| . 

The first reason can be characterized heuristically by saying that each eigenvector with a 
complex Uj is a discrete resonance, and there may be "background" under the resonances. A 
more formal statement is that there may be contributions to the Green's function other than 
from the discrete poles uj in the complex frequency plane; this possibility will be further 
explained below, but will not be the focus of the present paper. 

The second scenario is more intriguing — even though the discrete poles may give all the 
dynamics, the associated eigenvectors may nevertheless be incomplete. The simplest case 
is for two (in general any M > 2) discrete poles to coalesce into a double (in general Mth- 
order) pole, e.g., upon tuning some system parameter(s). We shall show below that, at least 
for 1-d wave systems, an Mth-order pole is not associated with M degenerate eigenvectors, 
but with only one. With the number of eigenvectors reduced by M — 1, must become 

incomplete, and H can only be reduced to block form (in this case an M x M block) rather 



than diagonal form. This is known as the Jordan block problem [13]. If it were to occur 



the above formulas (|1.3|), (|1.7|), and ( |1.8| ), derived assuming an eigenvector basis, would not 



hold. Such problems can be handled as the limiting case of several nearby simple poles, but 
since some eigenvectors are "lost", the limit is singular and needs to be treated with care. 

For a non-hermitian operator, in general there is no positivity condition to guarantee 
that (fi\fj) 7^ for its left and right eigenvectors {p\ and \fj) corresponding to the same 
eigenvalue u>j, even if both these vectors are unique up to a constant; in Section [FS| this 
problem did not occur only because was orthogonal to all other eigenvectors \fk) (k ^ j) 
of a set which was supposed complete. If (P\fj) should vanish for some j, then the above 
formulas would also fail. This is the orthonormalization problem. (So long as (P\fj) ^ 0, 
one can adopt a normalization convention that it is unity.) 



C. Outline of paper 



In most applications, it is implicitly assumed that the Jordan block problem does not 
occur, and the BB formalism has not been extended to handle such problems if and when 
they do. In part, this is because phenomenological NHHs do not exhibit these problems in 
a natural and convincing way, through which a possible extension of the formalism could be 
explored. 
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Recently it has been shown that waves in certain 1-d open systems provide an exact 
realization of BBs fl4j| . In these systems, dissipation occurs by the leakage of waves out 
of the system (through the imposition of the outgoing condition at the boundary). Under 
some simple conditions, the dynamics is completely controlled by the discrete poles in the 
frequency plane [|l|,[nj and an analogy with conservative systems can be developed ||17|| , 
such that the formalism can be cast in the language of BBs. In terms of these physical 
realizations, we shall show that the question at hand can be reduced to investigating the 
existence and properties of higher-order poles in the frequency-domain Green's function. 

In Section [TT], waves in open systems are described, focusing on the representation of the 
dynamics by the discrete poles [|18|] and the relation to a BS. We shall deal with both the 
wave equation and the Klein-Gordon equation. The duality transformation, central to the 
BS formalism, is shown to arise naturally from the dynamics. In Section ffTTl we construct 
examples where higher-order poles do occur; the examples correspond to the critical damp- 
ing of an oscillator. We then show that these higher-order poles relate to the Jordan block 
problem ("losing" one or more eigenvectors from a finite subspace). But since the dynamics 
is contained in the Green's function, a systematic study of the residue at a higher-order pole 
should reveal how the situation must be handled. We therefore first analyze, in Section [TV], 
the finite-dimensional subspace associated with the higher-order pole, i.e., the matrix rep- 
resentation of one Jordan block. Subsequently, in Section [V] we solve for the field evolution, 
which leads to an extension of the BB formalism with a generalized duality transformation. 
Then, in Section [VI], the perturbation theory associated with such a block is presented |19 



Unlike the degenerate perturbation theory of a conservative system, M — 1 of the basis 
vectors are not eigenvectors. Concluding remarks are given in Section |V1J| . 

The overall results of this paper are then (a) nontrivial examples of the Jordan block 
problem in some exact treatments of dissipative systems, (b) an extended BB formalism for 
dealing with such situations, and (c) a matrix representation of such Jordan blocks including 
their perturbation theory. 



II. WAVES IN OPEN SYSTEMS 



A. Wave equation 



We consider waves in 1 d described by 



p(x)d t - d x 



<f){x,t) = 



(2.1) 



on the half line [0, oo), with (p(x=0,t) = and <i>{x,t) approaching zero rapidly as x — > oo 
2H| . Let the system S be the interval [0, a], and the bath B be (a, oo), with p(x>a) = 1. 
Energy is exchanged between S and B only through the boundary x = a. The out going- wave 
condition d t (j)(x, t) = —d x (p(x, t) is imposed for x > a. 

The physical relevance of this model in describing strings 



21 1, electromagnetism [22 



and gravitational waves 
elsewhere |T|,[T5|ri|,[T3J2| . 

Everywhere below we restrict to p(x 
a density pT],p3|] or a dielectric constant 



as well as the mathematical formalism, have been discussed 



> 0, which is necessary both physically (p is 
and mathematically (the equation becomes 
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singular if p(xo) = for some xq). In particular, this means that p{x) can contain a 
(positive) ^-function, but not 5' or higher derivatives (for which equations such as ( p.3|) 
would be distributionally undefined anyway). Thus, has to be continuous, but <p' may 
have discontinuities. 

The eigenfunctions or quasinormal modes (QNMs) are factorized solutions 



(M) = fj(x)e 



-icjjt 



The spatial function satisfies 



[d 2 x + p(x)uj]}f J (x)=0, 



(2.2) 



(2.3) 



with fj(x) oc e luJjX for x > a. Introduce the conjugate momentum = p(x)d t <f), and the 
two-component vector 



(2.4) 



In terms of this, the dynamics can be cast into the Schrodinger equation idt\ 
the NHH 



H 



p{x)~ l 



2 







H\6) with 



(2.5) 



The identification = pd t (f) follows from the evolution equation ||25|| . In this two-component 
form, the eigenvectors are 



\fj) = 




U 

-iUjpfj 



(2.6) 



We shall define the duality transformation only after the dynamics has been discussed. 

The hamiltonian fl2.5|) is non-hermitian, as will be obvious once the appropriate inner 
product Q2.17] ) is defined. Nevertheless, on the "universe" [0, A] (with a nodal condition at 
x = A — > oo) a complete real spectrum of "universe modes" is guaranteed to exist. The 



latter are given by (2J3) for ±c<jj, where fj are the eigenfunctions with eigenvalue to 2 - of the 
operator — p(x)~ l d 2 , which is hermitian and positive on the "one-component" space with 
inner product (u,v) = J^dx p(x)u*(x)v(x). This indirect construction of the spectrum of 
( p.5| ) fails on the interval [0, a] (the energy of which is not conserved), however, since the 
very definition of the out going- wave condition <p(a + ) = —<f)'(a + ) is not possible with only 
one component. Consequently, the QNM frequencies will have a nonvanishing Imujj < 0. 



B. Green's function and poles 

The dynamics of these open wave systems is best discussed in terms of the Green's 
function G(x,y;t), viz., 



(f){x,t) 



G(x, y; t)4>(y) + d t G(x, y; t)p(y)<f>(y) dy 



(2.7) 
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where <fi and are the initial values of the wavefunction and conjugate momentum. The 
behavior of G is completely given by the contributions from the singularities in the frequency 
domain of its Fourier transform G(x,y;u). It can be shown [15] that under two conditions, 
viz., (a) p(x) has a discontinuity at x = a to provide a natural demarcation of the system 
from its surroundings (discontinuity condition), and (b) p(x) = 1 for x > a, so that outgoing 
waves are not scattered back into the system (no-tail condition), the only contributions will 
be from the isolated poles associated with the eigenvectors, and if all poles are simple then 
G can be represented as 



G(x,y;t) 



7 \fjlfj) 



(2.8) 



The normalization factor (fj, fj), which will turn out to be the same as (f-*\fj) and therefore 
central to the issue at hand, will be derived and discussed in detail below. Note that the 
numerator of G goes as fj(x)fj(y), not as, e.g., fj(x)f*(y) (which indeed would violate the 
symmetry of G in x and y). Thus (fj, fj) will have to go as fj, not |/?|; this makes a crucial 
difference, since QNM wavefunctions in general are not real. 
The derivation of the pole contributions is given in 



151, and the main elements will be 



reviewed below as preparation for Section One starts from the defining equation for the 
frequency-domain Green's function 



D(uj)G(x,y;u) 



d 2 x + p(x)uj 2 G(x,y;uj) 
-5(x - y) . 



The solution for 0<x<y<ais given explicitly by 

f(x,u)g(y,u) 



G(x,y;uj) 



W(u) 



(2.9) 



(2.10) 



where / and g are solutions to the homogeneous equation D(uS)f = D(u)g = 0, with / 
satisfying the left boundary condition f(x=0,uj) = 0, and g satisfying the right boundary 
condition g(x,u) oc e luJX for x > a. Their wronskian is W(uS) = f'g — fg', so that the 
combination ( |2.10| ) is independent of the normalization of / and g. 

One next writes an inverse Fourier integral to obtain G in the time domain, and distorts 
the integral from the real u axis to a large semicircle in the lower half-plane. Under the 
conditions stated, the contribution from the large semicircle vanishes and there are no cuts 



(or other singularities) associated with the tail of p [0,0 ; one is then left with the residues 
at the poles, namely the zeros of W ||26|| . (If the conditions do not hold, there will be non- 
pole contributions from the large semicircle and/or cuts, and the discrete eigenvectors will 
be incomplete in the first sense mentioned in Section |, leading to transients for short times 
and power-law tails for long times ||27|| . This scenario will not be pursued further here.) 

At a zero Uj of W, the functions / and g are linearly dependent: fj(x) = f(x,u>j) = 
Cjg(x, ujj). Thus fj satisfies both the left and right boundary conditions, and is an eigenfunc- 
tion. Assuming for the moment that Uj is a simple zero, the residue is related to dW(u)j)/dw, 
and it is straightforward to show that 



-a 



dW(uj 
duj 



2uj / p(x)fj{x) 
Jo 



dx + ifj(aY 



(2.H) 
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while in the numerator we have f(x,u)g(y,uj) = Cj fj(x)fj(y). The representation 
for the Green's function then follows trivially. In particular, the calculation shows why 
(fj, fj) oc fj instead of \fj\. The definition (|2.11|) of (fj, fj) is valid for eigenfunctions only, 
and will be shown to be a special case of (|2.13|) . 

This derivation and especially ( |2.11| ) make it clear that the orthonormalization problem 
is the same as the possibility of higher-order zeros in W, i.e., of non-trivial Jordan blocks, 
as was already apparent in the general analysis of Section ITB. 



C. Duality transformation 



While the duality transformation is simple and has been presented elsewhere ||14j| , it is 
useful to emphasize how it comes about naturally from the dynamics. For this purpose, we 
and find that \4>{t)) (cf. 



put (EJ) into 



is given by ( |1.8|) , with the coefficients [28 
fMky) + fiWiv)] dy + /i(a)0(a)) . (2.12) 



All reference to initial data on the "outside" x > a has been eliminated. Physically, the 
outgoing condition ensures that the outside data do not propagate in; mathematically, both 
the initial data and the Green's function are annihilated on the "outside" by dt + d x . 

The projection formula (|2.12|) suggests the definition of a bilinear map f29j between two 
vectors 



(ipX + 4>X) dx + tp(a)x(a) 



(2.13) 



It is seen that the definition ( |2.13j ) is consistent with (|2.11| ). Using this notation, we can 
now write the projection formula ( f2.12| ) compactly as 



(frJ)) ' 



(2.14) 



The map (|2.13|) has two important properties: (a) in the integral it cross-multiplies the 
two components, and (b) there is no complex conjugation of the first argument. These 
suggest the definition of a duality transformation T> = J 7 , with T a flip map that (a) 
interchanges the two components and (b) complex conjugates (cf., e.g., PD[), viz., 



T 



V>2 



r 2 



(2.15) 



in terms of which the bilinear map can be related to the standard inner product between 
two-component vectors 



{fax) = Wlx) 



(2.16) 



where the latter is defined as 



<Clx> = / (Cx + Cx)dx. 



(2.17) 
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In showing the equivalence ( [2.16| ), one has to use the outgoing property to collapse the 
integral over (a, oo) in ( 2.17|) to the surface term [17|. That is, ( p. 16 ) only holds when either 
\x) and or J-\x) an d |0, are both outgoing, and only in these cases we will use ( [2.17| ) 

in this paper. Thus, \x) and \() in general belong to different spaces, and hence the structure 
is not precisely that of a Hilbert space. Since this latter issue is unrelated to the possibility 
of having higher-order poles, however, these functional- analytic details will be dealt with 
separately fl3~T| . 

Two symbols T> and JF have been used in anticipation of their inequivalence in the case 
of higher-order poles. However, in the general case T remains defined by ( f2.15| ); cf. Q5.7| ). 
The property 



1731 ) is equivalent to the symmetry of H under the bilinear map: 



(2.18) 



To verify fl2.18| ), one needs to integrate by parts; the surface term generated exactly cancels 
the surface term in ( [2.13|) . The "orthogonality" of QNMs, 



(fj, fk 







(2.19) 



now follows in an immediate transcription of the standard proof. 

What remains then is to show that higher-order zeros in W can indeed exist (Section |Tl \ 
and to generalize the above formalism to those cases (Sections ^VJ-jVJ). 



D. Klein Gordon equation 



Also of interest in this context is the Klein-Gordon equation 

~d 2 t -dl + V{x)\ <f)(x,t) = 0. 



(2.20) 



Among other things, this equation describes the propagation of linearized gravitational waves 
on the curved background of a black hole p2 |. 

Essentially the entire formalism for the wave equation applies, mutatis mutandis |16 . 
First, the discontinuity condition now refers to V(x), and the no-tail condition requires 
V(x) = for x > a. Everywhere else we replace p(x) i— > 1, —d 2 x \— > —d\ + V(x). An example 
will be given below in terms of this equation. 



III. HIGHER-ORDER POLES 
A. General remarks 



For 1-d conservative systems, e.g., ( |2.3j ) with nodal boundary conditions at x = and 
x = a, W can only have simple zeros since in ( p. 11 ) the surface term is now absent, while 



the integral is positive definite up to an overall phase. Thus, the eigenfrequencies of conser- 
vative systems have a finite spacing Auj, and can be labeled by the number of nodes of the 
corresponding eigenf unctions. Hence, it is by no means obvious that higher-order poles can 
exist in the case of outgoing waves. For this reason, we will first demonstrate some examples 



S 



of higher-order poles (cf. p3[ ) before studying the extension of the BB formalism. For sim- 
plicity, here we shall concentrate on second-order poles. The demonstration of the existence 
of third-order poles and the discussion of other possibilities are deferred to Appendix A. 

Second-order poles can be obtained by allowing two first-order poles to merge. The con- 
verse is also true: any second-order pole, when suitably perturbed, splits into two first-order 
ones, as demonstrated explicitly in Section [V]]. One is thus permitted to think generically of 
the coalescence of poles, and this leads to the concept of the "loss" of eigenvectors. Consider 
two nearby first-order zeros of W, associated with two distinct eigenfunctions, say fj and 
fk- But when the poles merge, there is only one eigenfunction — in contrast to the case of 
degeneracies in conservative systems. To see this, it suffices to notice that for a given uj, the 
boundary conditions f(x=0,u) = and f'(x=0,uj) = 1, say, uniquely specify one function. 
The "lost" or "missing" eigenfunction makes the possibility of higher-order poles all the 
more interesting. 

It is best to look for second-order poles on the imaginary axis in the uj plane. First, apart 
from some overall factors of i, the problem is purely real and easy to handle. More physically, 
as the system parameter(s) are tuned, it is "unlikely" that two poles in the complex plane 
would collide — not only would this require the simultaneous tuning of two parameters, but 
one also expects level repulsion [3~3|. However, eigenvectors of the dissipative system ( |2 . 1| ) 
exist in pairs, with frequencies uj and — uj* lying on the same horizontal line. It would require 
the tuning of only one parameter to make them collide; when they collide, they must do 
so on the imaginary axis. In fact, we expect that after they collide, the two poles will 
move apart along the imaginary axis, in exact analogy to an oscillator going through critical 
damping. This scenario is exemplified in both models shown below, and it remains an open 
question whether higher-order poles can exist off the imaginary axis. 



B. Example in the wave equation 

With these remarks, we now look for a double zero of W in the case of the wave equation 
for uj = —i'j, with 7 positive. The differential equation (suppressing the mode index j |35) ) 
then becomes real: 

~d 2 x -p(x) 1 2 ]f(x) = (3.1) 

and the eigenvalue condition is /'// = 7 at x = a + , which ensures that W(— 27) = 0. For 
uj = —27 to be a double zero, we also need (/, /) cx dW/du to vanish: 

•(/, /) = 2 7 r + p(x)f(x) 2 dx - f(a) 2 = . (3.2) 

•> 

By using (|3.1| ) to express p(x)f(x) in terms of /" and then integrating by parts, this condition 
can be recast as 



«'(/./) = -- r f' 2 dx + f(a) 2 = 0. (3.3) 
7 Jo 



The last term has just been reversed on account of the surface term —2f (a) 2 produced in 
the integration by parts. 
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Interestingly, p(x) does not appear in ( |3.3|) , and this is central to the construction of 
examples, as follows, (a) Choose any function f(x) satisfying f(x=0) = 0, f'(x=0) > and 
f"(x) > 0. (b) Use ( |3.3| ) to determine 7. (c) Put these back into ( |3.1| ) to find p(x), which 
is guaranteed to be positive. 

There is however one further subtlety. Such a construction gives /'(a"), and also f'{a + ) = 
7/(0); the difference between these two must be attributed, through (|3.1|), to p(x) — ■ ■ ■ + 
pS(x — a), with p = 7~* — f'(a^)/('~f 2 f(a)). One must check that p > 0, i.e., that 



a 

12 



2 / f u dx > f(a)f'(a-) . (3.4) 
Jo 

This condition is nontrivial, and for instance violated for f(x) = x + ax n and for some a if 
n > 5. 

Yet examples abound, e.g., 

f(x < 1) = sinh(Kx) (3.5a) 

K 2 

7 = K cotanh K H — (3.5b) 

sinh K 

p{x) = %d{\ -x)+ 2 K [ 2 5(x - 1) + 6(x - 1) , (3.5c) 
7 Z r y z sinn ii 

for any K > 0. Note that in this example always p(0<a;<l) = K 2 /^ 2 < 1, the case in which 
there is a zero-mode (i.e., a single simple pole on the negative imaginary u axis) even if 
the 5-term in ( |3.5c ) for p is absent [1^,36]. Incidentally, this case without the 5-term thus 



already shows that the open string model Q2.1|) can exhibit features not found in damped 
harmonic oscillators, as will become even clearer in Appendix A. 



C. Example in the Klein Gordon equation 



We next give an example for the Klein-Gordon equation; this example may appear more 
natural in that the system (i.e., V(x)) is specified in advance and not obtained as an answer. 
Recall the Poschl-Teller potential |37j 



V(x) = V sech 2 x . 



(3.6) 



The Klein-Gordon equation ( |2.20|) with this potential is exactly soluble ]37| , with the eigen- 
values given by 
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±JV -\-i{j + l) 



if V > \ 
if ^0 < \ 



(3.7) 



j = 0, 1, 2, . . . . Each pair of poles merge at the parameter value Vq = \. Again, the merging 
is exactly like an oscillator going through critical damping. 

This example may however be regarded as slightly unsatisfactory in one way: the po- 
tential has no discontinuity but does have a tail (i.e., V(x) does not vanish outside a finite 
interval \x\ < a), so the set of discrete eigenvectors would not be complete even when all the 
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poles are first-order. Nevertheless, we may consider a minor alteration: if V(x) is truncated 
at |x| = L for some large L, one would expect the eigenvalues not to be much affected. 
Actually this is not the case in general P5[ , but here we shall not go into this subtlety, 
which is not related to the present issue. It is sufficient for the present purpose that at least 
the first pair of modes j = are not too much affected by truncation, in the sense that 
the L — > oo limit recovers the pole position for the untruncated V(x). This pair of poles at 
critical damping then demonstrates a double zero of If in a context where the eigenvectors 
are otherwise complete. We have verified numerically the existence of this double pole for 
L = 5, V = 0.252279109 . . itu = 0.511109 .... 



IV. JORDAN BLOCKS 

In this Section, we focus on a single Mth-order zero u>j of W and the subspace associated 
with it. As already mentioned in Section [111 A|, in this subspace there is only one eigenvec- 



tor, and presently M — 1 other basis vectors will be obtained. This means that H is not 
diagonalizable, and we consider this M x M block in H. In this and the following Section, 
the discussion will be confined to the wave equation (|2.1|). 



Using the definitions below ( 2.10|) , the (position- independent) wronskian can be written 
as W{oj) = f'(0,uj)g(0,uj). Differentiating with respect to lu, one now proves by induction 
that d™g(0, cu)^^. = for < n < M — 1. This means that, up to this order in u — Uj, the 
functions / and g satisfy the same boundary conditions and hence can be normalized to be 
equal, i.e., if we define 

M-x 

/(x,w)= £ f jjn (x){u - u,j) n + 0[{u, - , (4.1) 

n=0 

so that fj, n {x) = (V 7l 0^3'/( rr 3 a ')lw=wj5 then we also have 

M-X 

g{x,u) = £ f j>n (x)(u - co.r + 0[(u - u 3 ) M ] . (4.2) 

n=0 

Now define the time-dependent functions 

f. n ( x ,t)^^d:[f(x,u;)e-^% =UJj 

' ■ 1 iuitA 



rv. 



-d^[g(x,u)e 



n (—it) 171 

E hn-m(x) [ —^e-^ (4.3) 



m=0 m 



for < n < M — 1, where the last line follows using (|4.2| ). These are not only outgoing 
solutions of the wave equation (since g(x,u)e~ tuJt is such a solution for any u), but satisfy 
the nodal condition at the origin as well (since f(x,u>)e~' luJt has this property for any uj). 
The momenta associated with the fj tU are 

fj,n(x,t) = p(x)fa n (x,t) 

= -ip(x)[ujf jtn (x,t) + fc^fat)] , (4.4) 
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so that the action of the hamiltonian is H\fj tTl ) = u)j\fj >n ) + |/^„_i), with = 0. 

For fixed j, the functions fj, n ( x , t) (n = 0, . . . , M — 1) are linearly independent (in the 
sense that a non-trivial superposition cannot vanish identically in x and t), as is obvious 
by looking at the highest power of t in each (the coefficient of t n e~ %Uit in fj >n is oc fj(x), 
which by definition does not vanish identically). Therefore the initial data \fj, n ) have to be 



independent as well ^0|, otherwise one would have a vanishing superposition evolving into 
a nonvanishing function. 

Thus the set {\fj, n )} n =o' 1S a basis, in which the hamiltonian reads 





( w i 


1 


■ 


■• \ 









1 • 


■• 


H = 










. 




V o 





■ 





(4.5) 



Each entry in ( |4.5[ ) is really a 2 x 2 matrix (cf. ( |2.5| )) operating on a two-component vector 
such as (|2.4|) , giving H the structure of a 2M x 2M tensor product. However, the very fact 
that {\fj,n)}^=s is a basis for H means that each of these 2x2 matrices is oc 11. 

While it is guaranteed that H can be cast into the so-called Jordan-block form (fOj) 
in a subspace with precisely one eigenvector (e.g., [flip, we have now established a basis 
with respect to which this is indeed the case, and related this basis to the solutions f(x, u), 
g(x, uj). The basis is not unique, since a rescaling f(x, uj) \— > J\f(uj)f(x, to) mixes the / Jin , with 
only fj = /j.o remaining invariant up to a prefactor; in fact, this rescaling is readily checked 
to generate precisely those basis transformations which leave the form (|4.5|) for H invariant. 
A further — essentially unique — specification of the basis will be made in Section M. 

Vectors from different blocks are "orthogonal" under the bilinear map ( 2.13Q , i.e., 



(fj,m fk,r 







(4.6) 



The proof proceeds by induction with respect ton + m. The case n + m = is the standard 
one of eigenvectors given in ( |2.19| ) . Now consider 



^j\fj,ni fk,m) 4" \fj,n—li fk,m) \H fj,ni fk,m) 

^k\fj,m fk,m) 4" \fj,m fk,m—l) ■ 



(4.7) 



On both sides, the second terms vanish by the induction hypothesis, so one is left with 
(ujj - LO k ){fj,n, fk, m ) = 0, proving ( [OP . 



V. GENERALIZED DUALITY TRANSFORMATION 

Having obtained the extra non-eigenvector solutions fj jTl (x,t) (n > 1) in the previous 
Section, we now investigate how these enter into the field expansion. In doing so we shall 
consider all the poles simultaneously, so that the block size M acquires an index j. As 
with the Jordan-block form of H (cf. the remark below (|4.5|) ), a basis dual to {\fj, n )} n =o 
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is known to exist on general grounds (see Appendix B for details), but it remains to find 
its explicit form and, if possible, to choose the original basis (which still is subject to the 
freedom pointed out below ( |4.5| )) so that the ensuing expressions will be as simple as possible. 

Our starting point is (|2.7|) , with G(x,y;t) = f(du/2Tr)G(x,y,Lj)e~~ luJt , where G is given 
by ( p.lOp . The wronskian has an Mjth-order zero [pEf| 

W{u) = Wj, Mj {w - ouj) Mj M(uj) 

= W hMj {u - Uj )"> + 0[(u - uj)^ 1 } , (5.1) 

with M.{ujj) = 1. We now make use of the remaining freedom f(x,u) \— > J\f(u)f(x,u)), 
taking = AA{uj)^ 1 ^ 2 + 0\(uj — u>j) Mj ] (which is analytic in a neighborhood of u> = ujj) 

and similarly for g, preserving ( |4.2[ ). After this transformation we have [i2~ 



W(u) = W jiMj (co - Uj ) M * + 0[(u - uj) 2 "'] , (5.2) 



where we draw attention to the order of the error term. Eq. ( |5.2|) will greatly simplify the 
formulas below. The contour integral for G is now straightforward, leading to 

p -iujjt M j~ 1 n (_jf\Mj-l-n 

G(x,y,t) = E^r- E E f^vlf, u) \ nV (5.3) 

for t > 0, which we assume throughout this Section. By symmetry, (|5.3|) also holds for 
< y < x < a even though this is not the case for the original ( |2.10|) . If Mj = 1 for all j 
this agrees with the known result for simple poles, which follows by combining (|2.8|) and 
( P-llj ). In general, one can rewrite 

G(x,y;t) = ^2— J2 fj,M 3 ~i-n(y)fjA x ^) ( 5 - 4 ) 

j lW 3,Mj n=0 

in terms of the functions fj^ n (x,t) defined in (O). Insertion into ( |2.7| ) yields the time 
evolution 

1 M 3 -\ 

(f)( X , t) = -J2 77} E (fj,Mj-l-n , 4>) fjA X > t) ( 5 - 5 ) 

j W 0Mj n=0 

in terms of the bilinear map ( |2.13| ). In particular, this holds for <p = fk, m , in which case 
terms with j ^ k vanish by (|4.6|) , and the linear independence of the fk,n(t) discussed above 
(|4.5|) implies that the coefficients on both sides of (|5.5|) are equal, i.e., 



\fj,nt fk,m) ^jk^n+m,Mj — l 

W jMj . (5.6) 

Of course, this intra-block "orthogonality" relation is conditional on the normalization 

W^ n \ujj) = for Mj + 1 < n < 2Mj — 1, imposed in ( |5 . 2| ) ; without this normalization 
one has ( |5.13| ) below instead. In view of ( p.!6| ), the relation (|5.6|) leads one to define 

|^>) = V\f j>n ) = F\f j)Mi -i- n ) , (5.7) 
where T is the flip operation in (|2.15|) . Thus 
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(f j ' n \ff,n>) = -W h M 3 5 jr 5 nn , , (5.8) 
where the constant of proportionality Wj t M- is nonzero by definition (cf. (|5.1|) ), so that there 



is no normalization problem. Equation (5H) is a significant result; it shows that, unless 



Mj = 1 for all j, the duality map T> no longer coincides with T: T> changes the intra-block 
index n of \fj >n )- Since it is the flip map which obeys 

TE = R^T (5.9) 



note that the proof outlined below (|2.18 ) does not invoke any assumptions on the block 



structure of H), the relation (|1.3| ) in general is not satisfied by the operator V im- 
plicit in ( |1.6|) . Since it is an immediate consequence of (|5.9| ) that JF carries right into 
left eigenvectors and vice versa, the left eigenvector corresponding to \fj) = \fjfl) is 
(Ffjfil = (/ J,A/j_1 | 7^ (P ,0 \ f° r Afj > 1- While the left and right eigenvectors thus are 
orthogonal as stipulated in Section [TB|, this does not lead to orthonormalization problems 



as these vectors are not the dual of each other. 

Using ( |5.6| ) and ( |5.7j ), one may write the final result for the generalized bi-orthogonal 
expansion as 

G(x,y;t) =i±^ 2^ —r f ' ( 5 - 10 ) 

j n =0 \Jj,Mj-l-n 5 Jj,n) 

so that the time evolution reads 

W)) = T, E „Lf \ \fj,n(t))- (5-ii) 



The 1 1 limit of G(x, y; t) then yields the sum rule 



n=0 KJj,Mj—l-n ) Jj,n 



while in the same limit, ( 5.11]) is indeed seen to be of the form ( |1.6| ). 



Equation Q5.6Q is the generalization of (|2.11| ) and (|2.19| ), while our proof is a slight 



simplification even in the simple-pole case, cf. ||15|| . The representation (|5.10| ) generalizes 
and ( p.ll| ) extends ( |1.8|) and (|2.14|) . Also the simple-pole counterpart of ( |5.12| ) is 



already known [£L5[]. 

One may ask to what extent the basis we have obtained is unique. On the one hand we 
demand that in our basis the Hamiltonian have the block form (|4.5| ), and below this equation 
it has already been remarked that this forces the functions / 3jtl to be of the form ( |4.3| ) for 
some normalization of g(x,u)). On the other hand, to ensure the simplicity of formulas like 
( |5.11| ), we require that the dual to \fj, n ) be some .F|/j',n'}> where ( fO| ) then forces j' = j. 
Now for any normalization of / and g (i.e., temporarily abandoning ( ^4.2| ) and ( |5.2|) ) one has 



(fj,n, 9j,m) — _ Wj,n+m+l (5.13) 

as long as n, m < Mj — 1 (implying that the bilinear map vanishes if n + m < Mj — 2), 
as can be proved by operating with YH>=o[{ n ~ t)\(m + t + l)\]~ 1 d2~ e d^i +e+l \ UJ=UJ ' =U j j on the 
identity 
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cu 2 - to' 2 ) f pf(u)g(u') dx = [iLu'f(a, uo) - f(a + , u)]g(a, J) + /' (0, u)g(0, u') . (5.14) 



o 



Hence, the product (fj >n , fj,M -i-n) is always nonvanishing, and the only way of achieving 
bi-orthogonality by setting other products to zero as in (|5.6| ) is to normalize W as in Q5.2|) , 
which obviously fixes the \fj, n ) U P to one overall prefactor per Jordan block. Hence, our 
choice of basis is unique up to these prefactors. In fact, further analogy to the customary 
treatment of the simple-pole case results if one sets = — 2o>j for all j, so that 

(fj,n,fk,m) = ^j^jk^n+m,Mj-'i-- With this preferential normalization, which is convenient in 
applications jKJ, the freedom discussed in this paragraph is reduced further to one overall 
sign per Jordan block. 

In closing this account of the generalized duality transformation, let us return to the 
example of Section [111 B| with M — 2. Using either ( |4.1| ) or the integral representation ||40|| , 
one can find the "preferred" second basis function (i.e., the one for which (/o,i, /o,i) = p5| ) 
corresponding to f = f as in ( |3.5a| ) as 



ki(x) = f a (x) + f b (x) , (5.15a) 
K 

fa{x) = i — x cosh(Kx) , (5.15b) 

7 

f b (x) = -i (~ + tanh J fs:sinh( J fs:x) , (5.15c) 

where the conjugate momenta f a and /& are given by ([4.4|) . While (|5.15| ) is given here for 
reference and further use in Sec. 0, there does not seem to be a simple physical interpreta- 
tion of this result. The normalization occurring in the field expansion (|5.11 ) is evaluated to 
be 

K 3 

(f°\fo) = (/ 0,1 |/o,i) = (/o,/o,0 = — cotanhK. (5.16) 

7 



Since the contribution of /& to the product ( |5.16|) is oc (/o,/o), it vanishes. However, in a 



calculation in Section |V1 C| this term in / ,i wm be seen to be essential for arriving at the 
correct result. 



VI. JORDAN-BLOCK PERTURBATION THEORY 
A. Formalism for the generic case 

The use of BBs places dissipative systems into a framework very similar to that for 
conservative systems. As a result, it is straightforward to develop time-independent per- 
turbation theory, in effect by transcribing textbook results for conservative systems, which 
nevertheless apply to the complex eigenvalues and eigenvalue shifts |jT?|j4"4"|| . 

It will presently be investigated how this formalism, previously developed only for QNM 
spectra with simple poles, is modified if at least one eigenvalue Uj is associated with a 
nontrivial Jordan block. The ensuing splitting of the multiple pole into Mj distinct ones 
(in the generic case, to be defined below) is reminiscent of the lifting of a degeneracy by a 
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perturbation (typically breaking some symmetry) in hermitian systems; however, important 
differences need to be pointed out. In the first place, in the models here under discussion 
an Mjth-order pole is associated with Mj degrees of freedom, of which only one is an 
eigenvector, as has been emphasized before ||45|| . Secondly, in the Jordan-block case the 



splitting generically is governed by only one complex parameter, and as a result the Mj 
complex frequency shifts are not independent; in fact, all their relative magnitudes and 
phases are predetermined, and only the overall magnitude and phase depend on the details 
of the perturbation, namely on the one complex parameter; cf. ( |6.4j) . While this feature may 
seem unusual compared to conservative systems, it actually simplifies the calculation. 

While Jordan-block perturbations thus differ essentially from perturbations of degenerate 
levels in conservative systems, we will study the former by a method which can also be 
used for the latter: transferring part of the perturbing NHH H' to the unperturbed Hq 
and treating this part exactly, upon which the remainder of H' can be dealt with using 
conventional non-degenerate perturbation theory. 

Specifically, let H = H + XH', where H is assumed to have a known Jordan-block 
structure as described in Sections [TV] and |V], H' accounts for a change in density <5(p _1 ), and 
where |A| <C 1. For simplicity of notation, it is supposed that there is only one Mj > 2; the 
generalization to several Jordan blocks is immediate. Consider the splitting part of H'\ 

fff _ \f^-i)(f j ' M ^\ H> \mf j \ f(m 
s_ (/^- i i/i,M,-i) (fm • 1 ■ } 

This has only one nonvanishing matrix element, viz., 

(f^-i mfj) 

(6.2) 



a 



(fj,H'fj) 

{fj,Mj-li fj) 

cu]J a+ 5(p~ l )p 2 ffdx 

\fj,Mj-l-> fj) 



(6.3) 



We transfer this part to the unperturbed hamiltonian: H = if + Aifg, so that the remaining 
perturbation is H' = H'—H' s . The perturbation is said to be generic iff a ^ 0. This definition 
of a generic perturbation will be justified below, by showing that for sufficiently small A the 
matrix elements of H' effect only a higher-order correction compared to the splitting caused 
by a 7^ 0. First of all, however, it should be noted that the representation ( |6.3| ) for an 
infinitesimal 5p is evaluated to read a oc Jq 5pfJ dx. Using the method of variation of the 
constant (cf. the inner integrand in jE]]) and this latter representation, a ^ is seen to be 
equivalent to d\g(0,Uj) oc d\W(ujj) ^ 0. In other words, if a ^ the Mjth-order zero in 
the wronskian at u = ujj is split already in lowest order in A. 

The eigenvalue problem for Hq can be solved in the Mj x Mj block associated with ujj. 
For the characteristic polynomial in this block one has det(H — cull) = {ujj— u;) Mj ' — (— ) Mj ' Aa, 
yielding the eigenfrequencies as 

u j>n = ujj + se 27rni/M i (6.4) 
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n — 0,1, ... , Mj — 1), where s = M tfXa is an arbitrary but fixed choice of the root. Thus, 



3 

the splittings Cjj >n — ujj all have the same magnitude oc |s| oc A 1 /^ and are equiangular, 
i.e., their phases have constant differences. Both their magnitude and the overall phase are 
determined by a. The corresponding eigenvectors of H are 

Mi -i 

\kn)= E s m e 2 ~ /Mj U>>> ; (6-5) 

m=0 

since the higher-order pole has been split into first-order ones, their duals read simply 

\f' n )=m,n). (6.6) 

It remains to account for H', using conventional perturbation theory, by evaluating its 
matrix in the new basis. The validity of this procedure is not entirely trivial, since the 
transformation from the basis {\fj. n )} to {|/j, n )}> effected by the matrix P mn = s m e 2 ™ mi l M 3 ( 
is singular in the limit A — > 0. For a justification, denote the matrix elements of H' with 
respect to the old basis as H' nm = (P' n \H'\fj,m)/(P' n \fj,n), an d evaluate the inverse trans- 
formation (in fact a discrete Fourier inversion) as (P _1 ) m n = M~ x s~ n e~ 2 ' Knm% l M] . In the 
basis which diagonalizes H , the perturbation is then given as 



lp,n\fji\f. \ M 3~ l 

\J \ 11 \Jj,m/ V"^ fr>-l\ U' P 

\J \jj,n) k,£=0 



Ma-1 



-LeMern-nkWMjjj^i-k ? 



where it is crucial that the power s 1 M] does not occur since H' Mj _ l0 = on account 



of ( |6.1|) . Thus the matrix elements of XH' are 0{s Mi s 2 ~ Mj ), which means that the first- 
order frequency shifts due to H' are 0(s 2 ), small compared to the lowest-order splittings 
Aujj 7 n oc s. With energy denominators given by uj^ n — &j,n> oc s, higher-order shifts are 
smaller still by successive powers of s. Of course, matrix elements of H' connecting the 
block associated with the unperturbed ujj to other blocks, or connecting two other blocks, 
can be handled without difficulty. 



B. Nongeneric case 

If a = the leading behaviour is determined by other matrix elements, and the splitting 
of ujj can be partial or, depending on the scheme of calculation, occurs only in higher 
order ||19|| . We shall not investigate the general case, but instead give an example of a non- 
generic perturbation which is relevant to the discussion in Appendix A. Namely, if a 4 x 4 
Jordan block is perturbed by the operator 



H' 



/0 0\ 

a 



\0 a 0/ 



(6.8) 
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(in the basis {| /?>)}), the characteristic polynomial is found to be det(H — uA) = 
[(u — ujj) 2 — \a} 2 . Thus, the fourth-order pole is split into two second-order poles, and 
the latter do not undergo further splitting to any order. Of course, this treatment does not 
address the question whether for some Sp the perturbation H' can have the form (|6.8| ) for the 
open wave system ( |2.1| ), even if a fourth-order pole is assumed to exist. However, H' as in 
(|OD at least satisfies the fundamental symmetry (f j ' n \H'\fj^ m ) = (/ J '' M J~ 1- ' TO |if / |/ 3 - ) jv^_ 1 _ n ) 
(i.e., reflection symmetry with respect to the NE-SW diagonal within one block), which 
follows from ( |2.18| ) and (|T 



C. Example 



Returning to the example (|3.5|) of Section [III B| , and bearing in mind the NHH action on 



a two-component vector given by Q2.5| ) [jjfifl , we study the perturbation p i — > /? l (x) + 
X9(l-x). 



From its definition ( |6.2| ), and using ( |5.16 ) for the normalization, one obtains 



a = ^[KtanhK -smh 2 K] . (6.9) 

That is, if A > the frequency shifts ±y/a are purely imaginary, while for A < they are 
real and of opposite signs. In other words, as A is turned from positive values through zero 
to negative values, the poles move together horizontally in the complex plane, merge, and 
then move apart vertically, in accordance with the general observation made at the end of 
Section DTIX 



Proceeding to 0(A), conventional QNM perturbation theory gives the next contributions 
to the shift as (f > \H>\f ,o)/(f > \fa,o) and (j**\H'\fai) / (^{fai) respectively ||. By 
(gTp, both are evaluated as \{H' m + H' n ) + C(A 3/2 ) = ^ + O(A 3/2 ), where the last equality 
follows from the symmetry pointed out below ( |6.8|) . In the numerator of 



tji if 0,1, H / ) , . 

H oo = —r f tt~ > ( 6 - 10 ) 

1/0,1) Jo) 



the contribution of fb as in ( 5.15c|) (which does not contribute to the denominator in ( |6.10| ) 



cf. below (|5.16|) ) is seen to be oc a upon comparison with ( |6.2|) . For both split levels the 
next-order shift thus reads 

„, (/-g'/o) .(2 K 1 V ..„> 

= i- — tanhiT tanh 2 K 6.11 

7 V 4 6 cosh isT 4 12 J y J 

and is seen to be purely imaginary, so that if the double-pole zero-mode is split along the 
imaginary axis in lowest order (i.e., (9(A 1//2 )) the perturbed QNMs stay on this axis up to 
O(X), again consistent with the symmetry argument of Section |III A| . 

For a check, the perturbed QNMs can be also obtained directly from the wave equation 
Q2.3| ) together with the boundary conditions. For p~ 1 (x) having a constant value (p') _1 = 
•y 2 /K 2 + A on < x < a, the eigenvalue equation is found to be 
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in which one has to expand y/p 7 = K/j — (K 3 /2 / y 3 )X + 0(X 2 ) and u = — ry + ooiVX + 0J2X + 
CU3A 3 / 2 + C(A 2 ). In O(A ) and G{VX), ( gig ) is satisfied identically. In 0(A) one obtains 
(lui^/X) 2 = a, with a as in Q6.9Q ; in 0(A 3 / 2 ) one obtains ^A = Hq , with i?Q as in (|6.11D . 
Hence, there is complete agreement between the Jordan-block formalism of Section |V1 A 
and direct expansion of the wave equation. 

Finally, an example of a non-generic perturbation is furnished by changing K to K' — 
K + A in (|3.5cj ), namely by Sp(x) = d K p(x), where in the differentiation of p its implicit 
K- dependence through 7 as in ( fj.5b|) must also be taken into account. In lowest order, H' 
shifts the double pole corresponding to K to a double pole corresponding to K', and indeed 
Jo fipfo is found to vanish, as stipulated below (|6.3|) . Since beyond this leading order 
H (K) + XH' 7^ Hq(K'), the double pole will be split eventually. In line with the treament 
in Section VI A| , however, this is not pursued further here. 



VII. CONCLUSION 

A remark is in place on the relevance of the issue considered in this paper. The above 
and especially Section |VI] make clear that the set of systems for which non-trivial Jordan 
blocks occur is of measure zero in parameter space. However, this feature is shared with, 
among others, stationary points in the phase space of dynamical systems, critical points in 
phase diagrams (note the semantic coincidence with "critical damping" ) , and degeneracies in 
conservative quantum systems, all of which are worthy of study and are known to determine 
the global structure of a system's parameter space to a much greater extent than one would 
think at first sight. In the case of degenerate quantum levels, a further motivation is their 
relation to a system's physical symmetries. While the corresponding phenomenon does not 
yet show up on the level of this paper, further investigation reveals that two states can merge 
in the superpartner of a spatially symmetric open wave system p7 |. 



The existence questions raised in Section [III] and Appendix A are closely related to 
the problem of spectral inversion. In conservative systems, the classic inversion problem is 
to determine the system (e.g., p(x) for the wave equation or V(x) for the Klein-Gordon 
equation) given all the real eigenfrequencies Uj. The solution to this classic problem is 
well known [ 48[| . The analogous problem for open systems is to determine p(x) or V(x) 



from the complex eigenfrequencies Uj, or, more generally, from the singularity structure of 
G(x,y;uj) in the cu-plane. If, for example, G is specified to have poles of order Mj (say, 
Mj = 4) at Uj, does a corresponding p(x) or V(x) exist (at least for one in a class of such 
singularity configurations)? Assuming the general inversion problem for open systems (a 
topic for further investigation) to be tractable, at this stage each of the following scenarios 
seem conceivable. 

(a) The inversion algorithm indeed yields a p(x) with, say, a fourth-order pole or a pair 
of off-axis double poles in its spectrum. 

(b) The inversion problem turns out to have no solutions, yielding a non-trivial proof of 
the non-existence of these more exotic configurations. 
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(c) This particular set of singularities points to limitations in the inversion algorithm 
which otherwise might have been overlooked. 



Any of these possibilities would further the understanding of QNMs in open wave systems. 

Recently, for the case of simple poles we have carried out the second quantization of the 
open wave system ( |2.1| ) using QNMs ||56|| . The QNM expansion coefficients cij (cf. ( |2.14| )) 



emerge as the pertinent quantum degrees of freedom, in terms of which it is possible to 
eliminate the "outside" from the equations for the cavity evolution. This relevance to the 



quantum problem further motivates the study of the mode structure of (2T), and second 
quantization in the case for which this structure involves nontrivial Jordan blocks indeed 
turns out to be possible either by Hilbert-space methods or by exactly solving the associated 
path integral |43[ . 



In closing, it may be useful to place the present work into the following context. Many 
wave phenomena in nature can be described by an evolution equation idt\<p) — H\<f>) (see 
above ( |2.5| )) and a natural question is: what are the possible forms for H — leading to various 
types of time evolution for |0) — and how are these exemplified in physical systems? The 
most familiar examples are conservative systems, for which |0) is expandable in a complete 
basis of normal modes, in terms of which H would be diagonal with real eigenvalues. Our 
earlier work [ |H|[rT[ shows another realization: in a large class of outgoing wave systems, \<p) 
is again expandable in a complete basis (of QNMs), in terms of which H would again be 
diagonal, but with complex eigenvalues. The QNMs, however, are not orthogonal under the 
standard inner product; therefore it is convenient to introduce their duals as well, together 
with which they constitute a BB. The present work has identified and studied a further 
generalization pertaining to such open wave systems, for which H is not diagonalizable. 
In these circumstances we have shown that a well-defined Jordan-block structure emerges, 
involving a nontrivial duality transformation. 
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APPENDIX A: THIRD-ORDER POLES AND BEYOND 

As stated in Section |111 A|, at present it is not known whether higher-order poles can exist 



off the imaginary axis. The analogy to damped harmonic oscillators suggests that they can 
not, but this does not lead to a proof directly since the harmonic oscillator picture itself is 
contingent on the QNM spectrum being simple, with at most a double-pole zero-mode. In 
fact, already at the end of Section ID [15] it has been remarked that the analogy is imperfect, 
and also intuition which is mainly based on the conservative and WKB limits could bias one 
against more exotic possibilities. 

A possible strategy for looking for double poles off the imaginary axis is to first construct 
a fourth-order zero-mode, which upon applying a suitable non-generic perturbation could be 
split into a pair of such double poles, cf. Sec. [Vl|. Thus returning to the case uj = — ry, ( |5.13|) 
shows that for a fourth-order pole one needs (f , io,i) = (/o,i ; /o,i) = [j35fl . Expanding / 01 
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by ( p.l|) and Q3.3|) respectively, one arrives at two functional equations for / — / ]49 



by its integral representation |40], using ( f4.4|) for its momentum, and eliminating p and 7 

(Al) 



i7 3 Wo,3 



1 dx 

T^M 



dyf"(y)f(y) 



- I (f') 2 dx = 
Jo 



dxf(x)f(x) 



dy 
P(y) Jo 



-4 



1 da; 

THxj 



dyf(y)f(y) 



dzf'{z)f{z) 







i 2 



(A2) 



where without loss of generality we have chosen a = 1, and where ( |Al| ) alone implies a third- 
order pole. Solutions to these equations are to be sought among the functions / satisfying 
conditions (a) below ( j3.3| ) and the inequality (|3.4|). 

Up to now we have only investigated third-order poles. The strategy is to seek a function 
/ satisfying (|A1| ) and, once one is found, to construct the corresponding p(x) using (|3.1|). 
We have done so within the class of functions f(x) = x + ax n (n > 2) mentioned below 
( p.4| ) p0| . It is easy to see that Wq >2 < both for small and for large a: (a) for small a, 
f" « so that the second term in (|A1|) dominates; (b) for large a, one can neglect the linear 
term and hence find that 7 2 W / o,2 = —a 2 n 2 (4n — 3)/ (2n — l) 3 < 0. However, there exist n for 
which W ,2 can become positive for a in an interval (cti, a 2 ) — so that W 0:2 — at a,\ and at 
a 2 — for example {a x ,a 2 ) = (2.059,3.8209) for n = 5, and («i,a 2 ) = (1.063096,8.30908) for 
n = 6. In both cases, the inequality (|3.4|) is violated at «i but satisfied at a 2 , implying that 



third-order poles indeed do exist. Besides being a stepping-stone in the seach for fourth- 
order (and hence off-axis) poles, this result in itself already justifies the general (i.e., not 



limited to Mj < 2) setup in Sections |TVHVI 



APPENDIX B: CONSTRUCTING DUAL BASES 

In Section 0, one is faced with the problem of calculating the basis dual to {\fj, n }}n=o 
There is a standard result for finite-dimensional spaces |^TJ, which however applies only 



when the dual basis is constructed within the space spanned by the original one. Here we 
are concerned with the original basis in V = £[{\fj jn )} n A } (£ denotes the linear span), 
but with the dual basis in a different space W. (Guided by the simple-pole case, we expect 
W = TV with T as in ( |2.15| ), but this property will not be used below.) 



Therefore we are led to consider the following problem. Let V be an M-dimensional 
subspace of a Hilbert space with basis {t^}^^ 1 , and let W be another M-dimensional 
subspace of the same Hilbert space. Under what conditions will there be a dual basis {w 11 } 
in W, in the sense that (w m \v n ) = 8™7 We claim that the necessary and sufficient condition 



is 



wnv ± = {o}, (bi) 

where V 1 - is the orthogonal complement to V. (For example, if the whole Hilbert space is 
3-d, and if V is the x-y plane, then W must not contain the z-axis.) 

For a proof, let {w n } be any basis for W. Then the duality of {w n } and {v n } is equivalent 

to 
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(wq\v M -i) 





f w° \ 




( wo \ 


) 









(B2) 



The necessary and sufficient condition is that the metric matrix on the LHS is nonsingular. 
Singularity of this matrix would mean that a nontrivial linear superposition w of the w m 
(i.e., a nonzero vector in W) is perpendicular to all v n , i.e., that w G V . This simple 
calculation not only proves our assertion but, in any finite-dimensional space, also gives a 
constructive algorithm. 

Returning to the system of outgoing waves under discussion, the contour-integral calcula- 
tion of the main text now in effect shows that (|B1[) is indeed satisfied for Vj = C[{\fj,n)} n =o } 



and Wj = J-Vj, solves (|B2j ) for this case, and extends this bi-orthogonalisation to the whole 
Hilbert space by showing that the latter equals ®jVj, with ( |4.6|) holding between different 
blocks. 
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